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THE LIFT DISTRIBUTION OF WIN&S WITH END PLATES * 

By W. Mangier 



SUMMARY 



This report describes the lift distribution on wings 
with end plates for the case of minimum induced drag (in- 
duced downwash constant over the span). The moments on 
the end plates are also determined. 

It is found that moving an end plate of certain 
length up from the symmetrical position, is followed "by a 
slight increase of the total lift. As a marked increase 
in the moments of the end plate about its attachment also 
results, the symmetrical end plates are most advantageous 
if a higher circulation on the wing only is contemplated. 
But, if the end plates. are to serve for producing hori- 
zontal control forces, a pronounced unsymmetrical arrange- 
ment will be preferred. 

The magnitude of the ensuing side force for end plates 
of any length and position is illustrated in figure 13. 
The ratio of coef f icientrs of side force and of lift can be 
interpolated from figure 14. For the application of the 
side force, an empirical formula (22) is given. 

The writer wishes to express his appreciation to Miss 
I. Lotz, for the many suggestions in the preparation of 
the study. 

NOTATION 



b, wing span. 

h, height of end plate. 

k 0 , ratio of height of the upper part of the end plate to 
semi span . 



*"Die Auf triebsverteilung am Tragflugcl mit Endscheiben . 11 
Luftfahrtf orschung, vol. 14, no. 11, November 20, 1937, 
pp. 564-569. 
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' NOTATION (Cont. ) 



k u , ratio of height of the lower part of the end plate 
to h/2. 

k = ratio of height of upper to height of lower part 

k u of the end plate (k = 1 for symmetrical end 
plates). 

- = kfi_t-i^, ratio of height of end plate to span. 

b c - 

I, coordinate (referred to -| = l) along the span. 

II, coordinate (referred to | = 1 ) along the end plate. 

£ = | + iT], coordinates in a plane through the vortex sur- 
face shed "by the wing (referred to b/2 = l). 

T, circulation. 

A, lift. 

c n , lift coefficient. 

Wj_, induced drag. 

c^^ , drag coefficient. 

w i5 rate of downwasb induced by the vortices leaving the 
ving . 

2wj = v; rate of downwash induced by the vortices at great 
distance behind the wing. 

V, air speed. 

p, air density. 

P 2 

q = g V , dynamic pressure. 
F = b t, wing area. 

S Q (or S u ), side force on upper (or lower) part of the 
end olat e . 
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NOTATION (Cont.) 

F 0 (or F u ), area of upper (or lower) part of end plate. 
c_ (or c_ ), side force coefficient referred to F n (or 

" u »„>. 

M 0 (or M^), moment exerted by S 0 (or S u ) about the 
point of attachment of the plate. 

I. INTRODUCTION 



tTingg with end plates (fig. l) are used on sailplanes 
and lately, also, on the tails of airplanes. 

The object of the present report is to ascertain the 
relationship of the circulation distribution over the 
wing and of the lift to the height and position of the end 
plate. The side forces and moments on the end plates were 
also determined. 

The investigation hereinafter treats the case of min- 
imum induced drag, i.e., the case where the rate of down- 
wash Wj_ induced at the wing "by the shedding vortices is 
constant over the span, and where the shedding vortex sur- 
face moves downward at a constant rate 2w^ = w like a 
solid "body (references 1 and 2). 

The case of symmetrical end plates (k = l) was 
treated according to F. Nagel (reference 3). For small, 
one-sided end plates (k = a>). I also had access to the 
unpublished calculations by F. Keune. 



II. LIFT AND MOMENTS 



Assume a plane placed at right angles to the stream 
direction through the shedding vortex surface at a great 
distance "behind the wing. Provided that the downward ve- 
locity of the vortex surface (which has the value of w = 
2w^ if Wj_ is the downward velocity created hy the vor- 
tices at the wing) is small- that is, no rolling up of the 
vortex surface takes place - this section has the form in- 
dicated in figure 2. Our problem can then he solved if we 
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treat the plane problem of the flow around the section of 
figure 2, through the vortex layer. To make the problem 
steady, it may "be assumed that the section which is at 
rest is approached from "below at the constant speed w = 



is known, the circulation strength at a point % of the 
wing or point Tl of the end plate (fig. l), is equal to 
the jump suffered "by the potential when changing at point 
i, corresponding to point | or point 1 + iT| , corre- 
sponding to T| , respectively, from one side of the vortex 
layer to the other (figs. 1 and 2). The section "being a 
streamline (reference 4), \|/ is constant. 

Thus, we have: 



where $i ' and $2 are the potential values for point 
t or 1 + ill). 

Then the lift "becomes: 



If the complex potential 



£\j = $ + of this flow 




(1) 



+ x 




(2) 



- 1 



With this 



+ 1 




or, making T 



dimension] e ss with 



w 



2 



+ 1 




The induced angle of attack becomes: 
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which with 



IH * tan Bl = -i = - W - = -i = (3) 
V 2V A c a 



= /^iIidT (4) 



J = 

o " 2 
gives the drag coefficient 

w i a ttTd 2J a TTh 3 

In this form the factor k = ~~ gives at once the 

2 J 

departure of the drag coefficient from that of the ellip- 
tical wing. For equal lift and dynamic pressure a wing 
of span d with end plates has the same induced drag as 

a wing of span without end plates. 

For the determination of the side forces and moments 
at the ond plate, the sense of rotation of the hound and 
free vortices is given in figure 3. Seen from above, the 
hound vortices of the right-hand end plate turn in the up- 
per part in positive direction; in the lower part, in neg- 
ative direction; so that the air speed V causes a force 
perpendicular to V and to the end plate. In the upper 
part this force is inwardly, in the lower part outwardly, 
directed. The forces created hy the right and left end 
plates cancel each other, hut exert a moment ahout the 
point of attachment of the end plates. 

For the side force S Q , in the upper part of the end 
plate, we have: 



s 0 = p v | / r (n) an (e) 



whe re 



o 

k k __ 2h 

0 " 1 + k h 



17 e write : 
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o 




j = / isni d Ti (?) 



and find from equations (2) and (6) in conjunction with 
equation (4), the ratio of side force to lift at: 



= £o ( 8 ) 
A 2J 



and correspondingly for the lower part: 



£u = £u (9) 
A 2J 



* - 1 + k h 

The side force coefficient c g is introduced through the 
■relation 0 

So = c Sq f V 2 F 0 

where F 0 denotes the area of the upper part of the end 
plat e . 

Assuming the wing to "be rectangular (F = t> t) with 
elliptical end plates 

(f 0 = ^ k 0 h t; F u = | k u D t) 

we find: 



where 



and 



f^o = J o F = £ o 8 _1_ 



-a 



2J 



2J tt k 0 



(8a) 



and 
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J u F 



£u 8 JL_ 
2J 2J it k u 



(9a) 



The moment M 0 of the side force S 0 , is given "by 



p.v - 4 - 



r (T| ) T| dTl 



(10) 



TTith equations (3) and (5) the moment ahout the upper 
part of the end rclate becomes: 



^ - J 
A * 



2J 



(11) 



and ahout the lower part: 



if we r>ut : 



M. 



u 



'u 



A 5 2J 
2 



(12) 



and 



- k. 



^ Ti dfj 



Ml rj dTj 
w $ 



(13) 



u 



Finally, equations (10), (8), and (ll) yield for the lever 
arm at the upper part of the plate, the formula: 



and at the lower "oart: 



M„ M n A J * . . 

T) n = r _ x o__ = -P. = -0- (14) 

S 0 A | S o J o 



* 

J u 

n u = j~ (15) 

d u 
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III. TRANS FORMAT ION 

The flow potential C|U (£) in the (£, T| ) plane is 
determined "by so transforming the (|, T| ) plane through 
an analytical function £ = f(z) into the z plane (z = 
x + iy) that the right £ half plane (£ > 0) is trans- 
formed into the upper z-half plane (y > 0) at the saine 
time as the right-hand part of the section and the imagin- 
ary £-axis changes into the real z-axis (fig. 4). Then 
the flow potential Cp ( z ) and consequently, also, dp ( O 
can he readily estahlished. 

The desired transformation could be easily expressed 
with the Schwa,r z-Chri stof f el transformation formula, which 
would give an elliptical integral whose parameters - the 
points of the z-plane for the corners of the £-plane, 
would have to he determined through the geometrical quan- 
tities of the £-plane. But as this would entail the so- 
lution of a system of five transcendent equations with 
five unknown quantities, a different way, an approximate 
method, is preferahle. The desired transformation is 
"built up in three stages: 

1st stage: £ = + JZ~ X (16) 

The half plane I > 0, to which we can confine our- 
selves for reasons of symmetry, goes into the z '-plane 
(figs. 5 and 6). 

The straight pieces BC and ED "become the parabol- 
ic arcs B'C 1 and E'D'. It is permissible to assume 
that the ratio k of the distances BC and ED is ^ 1 : k= 
k 

— - - 1. In addition, since we treat only the cases where 

k Q < 1 (i.e., the upper end plate piece is less than half 

the span), the parabolic arc C'B'E'D 1 may be substituted 
by a circle. 

This circle is to be so chosen that its center lies on 
the real axis and passes through 3 1 and C'. Its equa- 
tion roads: 

x ,2 - 2ax ' + y » 2 = 1 - 2a 

whereby the center lies at the point a = - (l + i k 0 2 ) . 
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It also intersects the real axis in the point Gr » = L 1 = 
- (3 + k 0 s ) (fig. 6). 

ITow every point of the parahola is projected parallel 
to axis y* on the circle. (Any other projection would 
serve as well; we chose the most convenient one.) This 
means that the investigation tends toward end plates with 
small outward curvature rather than flat ones. This curva- 
ture, which depends on h/b, is negligibly small in our 
case . 

Point D« = $1 - k u s ) - 2i k u 

is replaced hy 

D» x = (1 - V 2 ) - 2i J 1 + \ (k 0 2 - k u 2 ) 



(3 + k c 2 ) z" 



2d stage: z » = (17) 

"0 



(4 + k 2 ) - ■*» 



A linear transformation transforms the point L* into 
point co • the circle therefore "becomes a straight line. 
Point H ' =:co becomes point H" = 4 + k 0 2 , and the end 

points C* And D' become the points 

c = i + i |a ( 3 + k 0 2 ) 

D- . 1 - i *sJ1±*l2 



The treatment thus reduces to the flow of a source 
and sink at point H" about the surfaces, illustrated in 
figure 7. Up to terms of higher order and a real factor, 
the potential in the vicinity of point H" = a 6 " has the 

form — - — » This is easily treated with the Schwarz- 

Christoffel formula. 

3d_st a ge: The outside of polygon A" B" C" D " E" P» 
Q.11 h " L" is transformed into the upper z-half plane 
(y > 0) (fig. 8). The effect of the approximate method 
is that the polygon now contains only two angles of magni- 
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tudc rr/2, and three angles of magnitude 2tt. The inte- 
gral following the application of the Schwar z-Chri stof f el 
formula therefore lends iteelf to evaluation without el- 
liptical integrals: 



z" = -• 



z 

1 P ( ?. - a 2 ) ( z ~_ a z)__(?_Z--llJ^l 

[i}i 'J J (z - a x ) (z - a 4 ) (z - a 5 ) 2 
a o 



a) x denotes a scale factor so chosen that point z" = 1 
corresponds to point a x . Oiving to the condition that the 
integral should "be real for real z > a 5 = 0 (fig. 8), 

and the further condition that the integral from aj to 
a 4 must he zero, it can even he expressed by an algebraic 
funct ion : 

■y + f-a _fj;g a 6 

ai a 4 z /a-, - z a 4 -z 

z" = 1 ± - — / -- 1 (18) 

+ a 2 a 3 a g y a x - a 0 a 4 - a 0 

a i a* a o 

Here : 

a Q < a x S a 2 % a 3 ^ a 4 < a 5 = 0 < a 6 =* - 1 

The prefix of the square root is to he so chosen that 
the upper half plane (y > 0) is transformed into the 
z "-plane; real z and z ^ ai are to carry the negative 
prefix, z ^ a 4 the positive sign, while for c.\ < z £r a 4 : 

■j_ + % a 3 a 6 

a i a 4 z / z - &x a 4 - z , . 

z" = 1 + i - — / - (18a) 

i + ^2_ a 3_ a 6 y a i - a O a ^ ~ a 0 

a l a 4 a O 

Five conditions exist for the parameters a 0 , &\ , a 2 , 
a 3 , a 4 (a 5 =0, a 6 = l). Two of these follow from the pre- 
viously cited demands that the transformal function for 
positive real z must te real, and that the length (C" D") 
must he equal to the sum of B" C" and D" B" : 
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a 2 + a 3 + a 6 = | (a x + a 4 ) 
a-i a 4 (a 2 a 3 + a 3 a 6 + a 6 a 2 ) = a 2 a 3 a 6 (a a + a 3 + a 6 ) 

From these equations a 1 and a 4 can he computed. 
A further deduction is that the function (18) in the point 
a 2 , a 3 , a 6 has the desired characteristics, .namely, that 
dz" 

— — — oecomes zero at these points. For this reason, there 

is no linear term at these points in the power series de- 
velopment for z"; (z" - a\") is proportional to (z- s.± 

(&l n is the imago point of a^ ; i = 2, 3, 6.) 

As the position of points C", D" , and H" is defined 
hy the position of C, D, and H in the £ -plane, the still 
free parameters a 2 , a 3 , and a 0 must he so chosen that 

point C" corresponds to point a 2 , point D" to point 

a 3 , and-point H" to point a 6 = 1. For a 2 and a 3 

(which alone occur in the first two conditions), wc can 
restrict ourselves to the range of 



as < 0; is = a 3 - " (1 + a s) - /l + a a + a 



3 

2 

The case of the wing without end plates (elliptical 
lift distribution ) is ohtained when putting 

a x = a 2 = a 3 = a 4 = - a 6 = - 1 (k 0 = k u = 0) 

while a 3 a 3 = a x a 4 = a 6 s = 1 gives the wings with sym- 
metrical end plates (k 0 = k u ) as computed "by F. Nagel 
(reference 3) hy a different method. 

IV. THE POTENTIAL 

Since our aim in the £ -plane at a great distance from 
the section is a flow parallel to the Tl-axis with a speed 
w (fig. 5), the development 

LpU) = w * (- i) i + .... (19) 
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is valid for the potential Cp ( i) in the vicinity of 

£ = CO* 

The transfornal function in the vicinity of £ = oo, 
that is, z = a 6 can he represented from equations (16), 
(17), and (18) in the following form: 

j = i __M — (1 + ... ) = i _JL_ (i + ...) (20) 

z - a 6 z - l 

Here 



and 



H = 



0 = 



with 



(3 + k n S ) (4 + k G 2 ) 



t§..-1..-± t&jUl. _I 

s , - a o a i - a o , . a a a 3 a s 

* 1 + 

a 0 a x a* 



(a 4 - a 1 ) a 2 a 3 



Oj. = -~ "a + - 

8 (a 6 - a x )" (a 6 - a 4 ) a x a 4 a f 



a 3 a 3 



a 1 a 4 



i±±2±±f- + _._2 a 6_-_Uj± a j2 ] (a =l) 



(a 6 -a x ) 2 (a 6 -a 4 ) 2 2a 6 (a 6 -a 1 ) (a 6 -a 4 ) - 



Then the potential "Cp ( z ) has, according to equations 
(19) and (20), the value: 

Cp( Z ) = w | = $ + iv|/ (21) 

Tor real ?. , = 0, i.e., a doublet with the real 
axis as a streamline. This result is anticipated (figs. 
5 to 8). No further terms occur in the development of 

cpu). 

Since ecuations (16), (17), and (18) afforded no sim- 
ple expression for z as function of £, several mathe- 
matical examples were computed. First the parameters a.± 
were determined. Then the correlated image points £ in 
the I- plane were computed for z values from the inter- 
val a 0 % z a 45 (point a 45 corresponds to point F ) 
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with the aid of equations (16), (17), and (18). The per- 
tinent potential Ep ( z ) = $(z) was obtained from equa- 
tion (20), after which the circulation r(£) and T (f\) , 
respectively, follows from equation (l). 

Mote.- The circulation on wings with plates (tail 
surfaces) of the form of figure 9 can he computed in ex- 
actly the same manner. The numerical results are to he 
published at a future date. 



V. NUMERICAL RESULTS 



Figure 10 illustrates various_circulation distribu- 
tions r(i) over the wing and r(T|) over the end plate 
as obtained with equation (l). It discloses that T(|) 
is substantially dependent on the height of the end plates 
i.e., on the ratio h/b, and that shifting the point of 
attachment of the end plate has no material effect on the 
total circulation. But near to the end plate itself, 
T({) is dependent upon the type of the plate attachment. 
According to equation (l), F is the difference of two 
potentials. They always proceed at the wing tip with ver- 
tical tangent, if the end plate does not protrude on the 
particular side of the wing; otherwise, the tangent is 
horizontal. In the first case the flow velocity from pros 
sure toward suction side is infinite at the plate attach- 
ment point; in the second case, it is zero. Thus, for an 
end plate at one tip (k Q = 0 or 1^ = 0 ) , T (| ) - has at 
point i = 1 a vertical tangent; otherwise, the tangent 
is horizontal. 

In figure 11 the value * = ~ (equation 5), that is 
the ratio — --- — is shown against k n for different 

a TTb 3 

values of h/b. The total circulation is smallest in the 
symmetrical case (k 0 = Jc^, k = 1 ) , (equation 4), but 
shifting of the point of attachment of the end plate re- 
sults in no appreciable change of ~ = \ (as anticipat- 
ed from figure 10), 



TT 



The point of application T] 0 of the side force S Q 
at the end plate as well as T| was determined from equa- 
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Tj was approxi- 
k u . 

(22.) 

The value 0.4 is an average value strictly applicable 
to symmetrical end plates (l: = l) only. But the discrep- 
ancy (0.39) for all other values, even k = » , is so 
small, that 0.40 applies to all k values. 

Figure 13 shows the S Q /A ratio against k 0 with 
k and h/b as parameters. The connection between S Q /A 
and k Q is the same as for S u /A and k u , if l/k is 
substituted for k in figure 13. In agreement with the 
fact that wings with end plates at one tip have a total 
circulation greater than wings with symmetrical end plates, 
the side force (per unit of length) with end plates at one 
tip exceeds that of those with symmetrical end plates.. 
Even the quotient S/A (for equal k 0 ) for end plates at 
one tip is greater than for symmetrical ones. 

The use of un symmetrical end platen therefore affords 
an increase in circulation on the wing (and consequently, 
a lift increase for a certain angle of attack) "but only at 
the expense of considerably greater side forces at the end 
plate. 

To illustrate: For J = 0.3, a shift of the end 

b 

plate from symmetrical to one-sided attachment, results in 
a 2-percent increase in circulation, but the side force on 
the upper part of the end plate increases almost three 
times, the arm at which it applies twice, resulting in a 
moment approximately 500 percent greater. 

To simplify the interpolation, we have plotted in fig- 
ure 14, the ratio c 0 /c_ against k n with k as param- 

oter. As with k 0 , both S Q and F 0 approach zero; the 

ratio c_ /c_ itself yields finite limiting values for 
o " 

k Q = 0. The relation "between c s /c a and k u is the same 
as between c s / c 0 and k Q , if k is replaced by l/k 
in figure 14. It will "be seen that c s '/c a can be approx- 



tions (14) and (15). It was found that 
matcly proportional to k 0 and T| u to 

k o k u 
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imately represented as linear function of k D for all 
values of k: c So / c a = A + B k Q . The A and B con- 
stants, being dependent on k, can he defined from fig- 
ure 14. 

According to the foregoing, it is advisable for 
structural reasons to employ symmetrical end plates. If, 

of course, the end plates are to serve as controls, then 

greater horizontal forces are desired. In this case k u 
would be chosen as small as possible. 

Regarding figures 11 and 13, it should be noted that 
only the shown points have been computed, after which the 
curves k = constant and h/b = constant, were obtained 
by interpolation. This method is well justified in view 
of the smoothness of the curves. 

A comparison of the theory with measurements by 0. 
Schronk (to be published in the near future), on a wing 
with end plate at one tip (k =»; k = 0.092), is seen, 
in figure 15. The recorded c a and c g values are plot- 
ted against angle of attack a; c s was obtained from 

pressure records (curve D) . It also shows the measured 

values for the moment coefficient of the side force 
M 

c mj; = — against ct, where F 0 is the area of the 

end plate and M 0 the moment of the side force S 0 about 
the end plate. 

The theoretical c s values (curve Th. ) were computed 
from the measured c a values by means of the relation 
c s/ c a = 0.446 from figure 14. The theoretical curve is 
in very close agreement with the experimental curve - at 
least, for positive a. The theoretical moment coefficient 
c mE of the side force can be obtained from equation (22): 

which, for the case in point, gives: 

c m™ = °« 4 x 0.446 c a (curve Th. ) 

These values are also in very good accord with the values 
c nE (curve D) obtained from pressure records. 
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ITo "better agreement between theory and experiment can 
be expected, particularly for c a = 0 , "because the rela- 
tion = 0 following from the theory of infinitely thin 
flat plate is not fulfilled for c a = 0. This is probably 
due to the fact that owing to the finite thickness of the 
profile and of the end plate and the profile camber, even 
at zero total lift, suction forces occur at certain points 
which occasion an inwardly directed force at the end plate 



Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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Fige. 1,2,3,4.5,6,7.8,9 
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Figure 1.- Wing with 

end plates. Figure 2.- Section through 

vortex surface Figure 3.- Bound and 



leaving the wing with 
end plates. 
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vorticee. 



Figure 5.- The £ plane. 
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Figure 4.- Conformal trans- 
formation of the 
£ on the Z plane. 
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Figure 6.- The 2' plane. 
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Figure 7.- The Z" plane, 



Figure 8.- The Z plane. 
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Figure 9.- Wing with 
. <-"' upper plate 8. 
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Figure Induced drag of wings 

with end plates. 
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Figs. 12,13,14,15 
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Figure 14.- Ratio c 80 /c e 

against height 
and position of end plate. 




Figure 15.- Theory compared 

with Schrenk's 
measurements. 
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